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A METHOD FOR THE CALCULATION OF EXTERNAL LIFT, MOMENT, AND PRESSURE DRAG 
OF SLENDER OPEN-NOSE BODIES OF REVOLUTION AT SUPERSONIC SPEEDS 

By Clinton E. Brown and Hermon M. Parker 


SUMMARY 

An approximate method is presented for the calculation of 
the external lift, moment, and pressure drag of slender open-nose 
bodies of revolution at supersonic speeds. The lift, moment, 
and pressure drag of a typical ram-jet body shape are calculated 
at Mach numbers of 1-45, 1.60, 1.76, and 8 TO; and the lift 
and moment results are compared with available experimental 
data. The agreement of the calculated lift and moment data 
with the experimental data is excellent. The pressure-drag 
comparison was not presented because of the uncertainty of 
the amount of skin-friction drag present in the experimental 
results. It was found that the lift coefficient definitely in- 
creased with increasing Mach number, whereas the moment 
coefficient taken about the midpoint of the body and the drag 
coefficient decreased with increasing Mach number. The 
manner in which the method may be applied to slender bodies 
of revolution with annular air inlets is shown. The excellent 
agreement of the calculated lift and moment results with experi- 
mental data indicates that the approximate method may be 
reliably used for obtaining the aerodynamic characteristics of 
slender bodies that are required for efficient supersonic flight. 

INTRODUCTION 

Current proposals for the design of aircraft capable of 
sustained flight at supersonic speeds and utilizing the ram jet 
as a method of propulsion have established the importance 
of knowing the aerodynamic characteristics of slender open- 
nose bodies of revolution at speeds greater than the speed of 
sound. The lack of theoretical treatments and experimental 
data emphasizes the need for theoretical investigation of this 
problem to serve as a guide for future work and as a check on 
the reasonableness of current and future experimental results. 

The small-perturbation approximation was used in refer- 
ence 1 to deduce the wave drag and in reference 2 to obtain 
the lift and moment of slender pointed-nose bodies of revolu- 
tion. No fundamental analysis is known to have been made, 
however, of the characteristics of a slender open-nose body 
shape, such as that required by ram-jet propelled craft. The 
peculiarity of the problem, from general considerations of 
similarity, is that the flow pattern is two dimensional at the 
lip of the nose and approaches the three-dimensional pattern 
farther along the body. The present work extends the 


method of references 1 and 2 to apply to these slender open- 
nose bodies of revolution with supersonic flow into the nose. 
The result is a fairly simple method of numerical integration 
of the differential equation of the flow. As an illustration, 
the pressure distribution, wave drag, lift, and moment are 
calculated at Mach numbers of 1.45, 1.60, 1.75, and 3.00 for 
a typical ram-jet airplane body shape; and the lift and 
moment results are compared with the experimental data. 
It should be pointed out that the accuracy of the method, 
which assumes potential supersonic flow throughout the field 
and also assumes small disturbances, depends on the surface 
angles of the body and the Mach number. The error in- 
creases with either increasing Mach number or increasing 
surface angles. 

SYMBOLS 


x, r, e 

cylindrical coordinates 

X 

distance along r-axis measured from nose of 
body 

l 

length of body 

R 

radius of body 

8 

I 

Mach angle ^sin -1 

<t> (x, r, 6) 

perturbation potential 

4> i (r, r) 

perturbation potential for axial flow 

<k (?, ri 6 ) 

perturbation potential for cross flow 

r* 

axial velocity increment ( ^ ) 

Vr 

radial velocity increment 

V 

velocity in undisturbed stream 

a 

velocity of sound in undisturbed stream 

M 

Mach number in undisturbed stream {V/a) 

p 

density in undisturbed stream 

Ap 

incremental surface pressure due to angle of 
attack 

Pi 

local pressure 

P 

pressure in undisturbed stream 

7 

ratio of specific heats of air (1.4) 

a 

angle of attack, radians (except where other- 
wise noted) 

5 

angle between surface of body and r-axis 
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c t 

0 D 

c n 

u 


lift coefficient ^Lift/^ WR#*^ 
drag coefficient ^Drag j~V i TrR N 2 '^ 
moment coefficient (^Moment 


variable of integration 


£=a -,—Br cosh u 

jt—Xt—BRi 

4, =/(*), 

N 
n 


deg 


SUBSCRIPTS AND SUPERSCRIPTS 

refers to nose 

refers to nth integration station, summation 
variable 

refers to ith integration station, summation 
variable 
in degrees 

MATHEMATICAL ANALYSIS 

POINTED BODIES OF REVOLUTION 


The linerarized equation of motion of a nonviscous com- 
pressible fluid may be written for a cylindrical coordinate 
system: 




a) 


where <£ is the potential function assumed to represent the 
effect of a small disturbance set up by the slender bodies 
being considered. The problem is to find a solution' of 
equation (1) that will satisfy the known boundary condi- 
tions at the surface of the body. A general solution of the 
differential equation (1), when 1, for diverging waves 
has been found by Lamb (reference 3) to be, with a slight 
change in notation, 


where 


and 


where 


<£=S Qj* C03 sd-hP/* sin s6 

S 


^‘“(rdr)^ 0 



/»09 

<2o= J f(x—Br cosh u)du 

n co 

P 0 = J g(x—Br cosh u)du 


B=-JW=i 


( 2 ) 


The part of Lamb’s general solution corresponding to con- 
verging waves does not apply to the present problem because 
all disturbances originate on the body and diverge into the 


flow field investigated. Von Kdrm6n and Moore have 
investigated the problem of the resistance of projectiles and 
cones (reference 1) and have found a solution for the case 
of axial symmetry 


<f>= — f(x— Br cosh u)du (3) 

that can be seen to be a special case of the general solution 
with s= 0. In their analysis the body, in this case a sharp- 
nose projectile, was represented by a distribution of sources 
along the x-axis starting at x=0, the nose of the body. By 
a numerical method of integration, it became possible to 
write the equations for the velocity increments v r and o x , 

Vr n =-Bg (4) 


where 
and 

A ( =f(0t 

with the boundary conditions 


Pr B =^ Attcosk-KTY*) — cosir ‘(TV!*)] 


(5) 


rp ft X * _2i 


j,=x t —Bli t 


v T dr 

Y+v x ~£c ‘ ( 6 ) 

These three equations in three unknowns (A n , t' r# , v x J were 
solved at each station on the body for v r and v x . The pres- 
sures were then found from the Bernoulli equation in the 
form: 



Ferrari (reference 4) and Tsien (reference 2) have inde- 
pendently found solutions for the case of pointed bodies of 
revolution at small angles of attack. Their solutions showed 
that the potential could be expressed in two terms: the first, 
from equation (3), 

<f>i= — /i (x—Br cosh u)du 

is the solution for the pure axial flow already described ; and 
the second 

<fa=—B cos jt(x—Br cosh u) cosh u du (8) 

represents the cross-flow potential of an arbitrary distribu- 
tion of doublets along the axis of the body starting at tho 
nose of.the cone or projectile. The form of equation (8) is 
that the cross flow is from the direction 9=0, as shown in 
figure 1. 
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By neglecting the small effect of the axial flow on the lifting pressures, Tsien obtained for the pointed projectile 
of arbitrary shape the equations: 





The values K t in these equations are assumed to be constants for each interval of the step-by-step process. The moment 
coefficient of equation (10) is assumed positive for nosing-up moments, these moments being taken about the nose. 


OPEN-NOSE BODIES 


The flow conditions over an open-nose body differ from 
those of pointed bodies in that, for finite angles of the nose 
lip, the flow is two dimensional at the lip. This problem 
was not considered in references 1, 2, and 4 and the general 
solution should therefore be examined to determine its 
applicability to this special case. Lamb has shown (refer- 
ence 3) that a sufficient requirement for the existence of the 
general solution to the differential equation of motion 
(reference 1) is that/(x — Br cosh u) be zero for all values of 
the argument less than some arbitrary limiting value. The 
determination of f(x—Br cosh u) such that the boundary 
conditions at the open-nose body are satisfied assures the 
fulfillment of this general requirement. For the usual case 
of supersonic flow into the nose, the boundary condition re- 
quires the surface of the body to be a continuation of a 
cylindrical stream surface of radius R H in the undisturbed 
flow ahead of the body as shown in figure 2. The perturba- 
tion potentials, equations (3) and (8), therefore must be zero 
at the cylindrical stream surface ahead of the body. Sub- 
stituting ^—x —Br cosh u in equations (3) and (8) gives 


•=-jr 


’*- ar um 

-» VCr-^’-BV 2 


A cosfl f'-* 

r J-» V( x-Q a -B¥ 


It is obvious that the boundary conditions are satisfied by 
letting y 1 ($)=/ 2 (f)=0 for all values of £<2r 0 — BR 0) where 



the point (xo, R a ) is at the lip of the open-nose body. It 
then remains to dete rmine (£) and/*(£) for £)>ar 0 — BR 0 so 
that the body surface is a continuation of this stream surface. 
From physical considerations, _/)(£) and /j(f) may be re- 
garded as an axial distribution of sources and doublets, 
respectively, where £ is measured along the r-axis. Because 
the effect of a source or doublet can be felt only along or 
behind its Mach cone, the source distribution must begin a 
distance BR 0 ahead of the nose. This point is chosen for 
the origin of the coordinate system. (See fig. 1.) It must 
be emphasized that the source and doublet distribution 
determined by satisfying the boundary conditions at the 
stream and body surfaces shown in figure 2 does not represent 
correctly the flow inside that stream surface. This result 
corresponds to the physical fact that the actual supersonic 
flow into the nose does not affect the flow external to the 
body. The basic assumptions of potential flow and small 
disturbances are valid provided the slope of the body surface 
is small. Actually, for finite angles of the nose lip, a non- 
conieal shock wave is formed that causes a loss in total 
pressure and produces rotation in the field. 

Numerical integration of equation (12), with constant 
values of fi(£) assumed over the integration intervals, 
results in the same expressions for r r „ and v Xn as those ob- 
tained by Yon K&rmdn and Moore (equations (4) and (5)) 
These constant values of/i'(f) are determined by satisfying 
the boundary condition: 


Vr 

V+v z 


=tan5„ 


(14) 


where tan d„ is the slope of surface of the body at the nth 
interval of integration. 
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By following the method of reference 2, the lift and moment coefficients for small angles of attack based on the area of 

the nose may be written 4 r-„ ci+brh ^ 

C L = ^-cos dddRdx (15) 


SL ( x ~^~ BR *)^t cos8deRdx 


where l is the length of the body, R x is the nose radius, and the moments are taken about the midpoint of the hotly. 
By substituting the expression for bfafbx in equations (15) and (16), C L and 0„ become 


0 


The distribution function /j(£) must be determined by the boundary condition; thus, 


Va cos 0=(|jr) (19) 

for which the radial velocity is assumed to be normal to the surface. -A more rigorous boundary condition taking into 
account the slope of the body was given by Ferrari (reference 4). For small surface angles, however, equation (19) is 
within the accuracy of the small-perturbation assumptions. The expression 

Vbr ) r . R ~W Jo ^J(x-)y-B i R 2 (20) 

is integrated numerically for constant values of over the ith interval of integration to obtain the sum 

n 

(b-?) V ^[cosh- 1 ^.,") -cosh" 1 m + (T m *) V(7 , < _ 1 »)*- 1 - Tf S(TW~1 ] (2I) 

Substituting this equation in equation (19) gives 

n 

[cosh- 1 (T,-,-) -cosh ~ 1 (T t n ) + (T { -i n ) VdV^-T- T? VOY^l ] (22) 

B?K 

With the values of ‘ determined, equations (17) and (18) become 
s\ a 

2 >1 1 s ^ 1 V vcr.T.->j*=i — } t23) 








l+2BR N x^+x. 


n 

=*)<*,-*_,) «,+«.-.) § j C Vc^t-iT-i- ■ 


+ [Vf^-i"- 1 ) 1 -!- vcr?= 


In equations (23) and (24) the pressure used for a given 
integration interval is the average of the pressures at the 
beginning and at the end of the interval. This scheme of 
using average lifting pressure is particularly necessary in 
regions where the pressure is rapidly changing. The method 
does not give the pressure at the beginning of the first 
integration interval, that is, at the point n=0. It can be 
shown that, as the first interval approaches zero, the pressure 
at the lip (n=0) is obtained by letting the expression in 
equations (23) and (24) 


n—l 



have the value 0.5 when n= 1. 

METHOD OF CALCULATION 

If calculations are to be made for an open-nose body, the 
total number of integration stations chosen is a compromise 
between the amount of labor involved and the accuracy 
desired, due consideration being given to the limitation on 
the accuracy imposed by the basic assumptions of small 


disturbances and potential flow. In general, where tho 
pressures are changing most rapidly tho integration stations 
should be the most dense. 


The integration stations must first be chosen. (See fig. 3.) 



Fiquhb 3— Diagrammatic sketch to Illustrate Integration process for body selected 
for calculations. 
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The calculation for the lift and moment then proceeds as 
follows: 

The boundary condition, equation (22), is applied to the 
point n= 1, the summation reducing to a single term, where 


and 


T< 


, gi-ji 


TV i 


~ BR, 


Equation (22) then gives a value of B 2 Kj2Ya. Next, 
equation (22) is applied to point 2, n=2, containing now two 
terms. Substitution of the expressions 

ATT* 2i 

12 ~^R7 

rp 2 *^2 3 1 

1 1 — br 2 

rn i J-2 0 

J ° ~~BR t 

and use of the value obtained for B 2 Ki/2Va permits the 
calculation of the value of B i K x (2Va. 

By continuing the process, at each successive station one 
more term occurs in equation (22) involving one new K u 
which is then determined. With the values of B 2 K t /2Va 
determined, equations (23) and (24) are used to evaluate 
the lift and moment coefficients. It will be noticed that, 
when equations (23) and (24) are evaluated for n=l, the 
expression occurring after the last summation symbol is 
unobtainable. For reasons previously stated, the expression 
must be given the value 0.5. 

The procedure for calculating the drag pressures is similar 
to but somewhat simpler than that for the lift and moment. 
Equations (4) and (5) are evaluated at the point 77=1, the 
sums reducing to one term. By applying the boundary 
conditions (equation (6)) and using the known slope of the 
body dr/dx, the constant A x occurring in equations (4) and (5) 
is determined. Substitution of A t back into equations (4) 
and (5) gives the increment velocities at n = l. Bernoulli’s 
expression (equation (7)) then gives the pressure ratio at 
n=l. It is to be noted, since V=Ma, that actually values 
of Ai/a, v ,,/a, and v x Ja are determined. Proceeding to the 
point 71=2, one new A is involved in equations (4) and (5) 
that is determined by the boundary condition in equation (6). 
In the same way, the velocities and the pressure ratio 
are calculated at 7i=2. The process is continued at success- 
ive integration stations over the whole body. The calcu- 
lated values of the pressure ratio allow a pressure-distribution 
curve to be drawn and the drag coefficient to be evaluated in 
the usual manner. 


DISCUSSION OF RESULTS 

RESULTS OF CALCULATIONS 

Calculations were made in order to obtain the pressure 
drag, lift, and moment of a typical open-nose body. A 
sketch showing the dimensions, the integration stations, 
and intervals is given in figure 2. Calculated pressure dis- 
tributions at zero angle of attack for the Mach numbers 
1.45, 1.60, 1.75, and 3.00 are presented in figure 4. The 
pressure rise at the nose lip is approximately that which 


would be obtained over a two-dimensional wedge of the 
same angle. Ackeret's theory for small disturbances (refer- 
ence 5) gives for the pressure rise of a 3° wedge ^=1.147 at 

.3/= 1.45, whereas the pressure ratio on the lip of the open- 
nose body is about 1.140. This agreement is considered a 
reasonable cheek inasmuch as the pressure on the lip must 
be an extrapolation of the pressure distribution from the 
first point of the integration process, which must be a finite 
distance back of the lip edge. It can be shown that, as the 
size of the first interval approaches zero, the pressure at the 
lip becomes that given by the two-dimensional theory. The 
effect over the nose section of the size of integration inter- 
vals is illustrated by figure 5. The 17-point method was 
chosen because a greater number of points resulted in only 
a small increase in accuracy and a large increase in the labor 
involved. As would be expected, for the case of the straight 
conical nose, the pressure falls off from the leading edge and 
approaches the pressure of a cone of the same surface angle 
(reference 1). At the corners of the body, the pressure falls 
approximately in accordance with the Prandtl-Meyer 
relation (reference 6) for supersonic flow around a two- 
dimensional corner. On the center and tail sections there are 
positive pressure gradients that for actual flight conditions 
would tend to cause separation. The source-distribution 
functions corresponding to the pressure distributions of 
figure 4 have been plotted in figure 6. 

The incremental surface pressures giving rise to the lift 
and moment are shown in figure 7 for the Mach numbers 
1.60 and 3.00. At the higher Mach numbers the pressures 
decrease less rapidly over the nose section. The doublet- 
distribution functions at the four Mach numbers are pre- 
sented in figure 8. The curves of figure 9 show the theoretical 
variation of lift-curve slope, moment-curve slope, and drag 
coefficient with Mach number. The rather interesting result 
obtained is that the lift coefficient increases with Mach 
number. Experimental results presented in reference 7 
show an increasing lift-curve slope with Mach number for a 
pointed projectile. The fact that the doublet distribution 
for the open-nose body is similar to that of a typical proj’ectile 
(reference 2) leads to the expectation that the Mach n umbei 
characteristics of the two shapes would be similar. The 
center of pressure at the lower Mach numbers is ahead of 
the nose and moves back with increasing Mach number. 
The drag coefficient can be seen to decrease with Mach 

number but to a lesser degree than the „ law of two- 


dimensional wing profiles. 


VM 2 -1 


COMPARISON WITH EXPERIMENT 

A comparison of the calculated lift and moment coefficients 
with some experimental data obtained in the Langley 9-inch 
supersonic tunnel is presented in figures 10(a), 10(b), and 
10(c). The contribution to the lift of the internal air can 
be shown to be: 

ACi=2a 

This increment will appear at the nose and will therefore 
give rise to a moment: 

A C~=a 
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o Body with ferrule 
[experimental) 
a Body alone 
[experimental 

— Ferrule tail 

— Regular tail 




(b) M=1.60. 

(c> Af-1.75. 

Figues 10. — Comparison of calculated lift and moment coefficients with experimental data. 
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These increments have been added to the calculated results 
for comparison with the experimental data that include the 
effect of the internal air. The experimental lift and moment 
results did not go through the origin because of deviations in 
the stream flow and perhaps small, unaccounted for, tare 
forces. It is obvious that for bodies of revolution the data 
should give symmetrical curves about the origin. The 
experimental results were therefore made to go through the 
origin by subtracting from each point the value of the lift 
coefficient or moment coefficient at zero angle of attack. 
All coefficients are based on the nose area. 

For M=1.45 two calculated curves corresponding to two 
body shapes tested are shown. The two shapes, which are 
designated for convenience regular tail and ferrule tail, 
differ only in that a stabilizer-supporting ferrule was placed 
over part of the tail section. The two configurations are 
shown in figure 10(a). No difference can be observed in 
the experimental results for the two shapes, whereas the 
calculations show a small difference. A possible conclu- 
sion is that thickening of the boundary layer ahead of the 
shock at the trailing edge tends to make the effective shape 
of the regular body more nearly like that of the body with 
the tail ferrule. The better agreement of the ferrule-tail 
calculation is further evidence that this thickening of the 
boundary layer actually occurs. In any case, the calcu- 
lated results are in excellent agreement with the experi- 
mental data. At M= 1.60 and M— 1.75, only the regular 
shape was calculated, and, although the agreement with 
the experimental data is good, ferrule- tail calculations 
would probably have given even better agreement. The 
discrepancy at M—1.75 between the theoretical and experi- 
mental points near ±3° is probably due to separation 
phenomena. The agreement is good at small angles of 
attack. Because of the uncertain value of the skin-friction 
drag present in the experimental data, a drag comparison 
is not presented. 

ANOTHER APPLICATION 

The success of the method when it is applied to open-nose 
bodies suggests its use in the calculation of external lift, 
moment, and pressure drag of bodies of revolution having 
annular air inlets (fig. 11). The characteristics of a body 
of this type would be calculated as follows: The source. and 
doublet distribution from 0 to B. would be numerically 



determined by the methods of references 1 and 2, which 
would allow the calculation of the lift, moment, and drag 
contributions of the portion of the body from 0 to C. In 
order to calculate the lift, moment, and drag of the remain- 
ing external shell, the source and doublet distribution from 
0 to A:would be retained and a new source and doublet dis- 
tribution beginning at A and satisfying the boundary con- 
dition at the outer shell would be determined. In (his way, 
the influence of the forward portion of the body on the flow 
along the outer shell would bo fully taken into account. 
The flow at the annular inlet could bo fully determined 
from the original source and doublet distribution from 0 to 
B , which would allow the determination of the lift due to 
the internal flow. The moment would be closely approxi- 
mated by assuming this incremental lift to act at the lip of 
the inlet. 

CONCLUSIONS 

An approximate method was presented for the calculation 
of the external lift, moment, and pressure drag of slender 
open-nose bodies of revohition at supersonic speeds. The 
pressure-drag comparison' was not presented because of the 
uncertainty of the amount of skin-friction drag present in 
experimental results. The calculated lift and moment 
results at Mach numbers of 1.45, 1.60, 1.75, and 3.00 showed 
excellent agreement with the available experimental data. 
The excellent agreement of the calculated lift and moment 
results with experimental data indicated that the approxi- 
mate method may be reliably used for obtaining the aero- 
dynamic characteristics of slender bodies that are required 
for efficient supersonic flight. 

Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 29, 19/$. 
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